o 






a^ 



^ 



in 

o 



S 



Flow pattern transition accompanied with sudden growth of flow 
resistance in two-dimensional curvilinear viscoelastic flows 

Hiroki Yatou 

Department of Physics, Graduate School of Science, 

Kyoto University, Kyoto 606-8502, Japan 

May 11, 2010 

Abstract 



We find three types of steady solutions and remarkable flow pattern transitions between them 
in a two-dimensional wavy-walled channel for low to moderate Reynolds (Re) and Weissenberg 
(Wi) numbers using direct numerical simulations with spectral element method. The solutions are 
called " convective" , "transition", and "elastic" in ascending order of Wi. In the convective region 
in the Re-Wi parameter space, the convective effect and the pressure gradient balance on average. 
As Wi increases, the elastic effect becomes suddenly comparable and the first transition sets in. 
Through the transition, a separation vortex disappears and a jet fiow induced close to the wall by 
^0 ■ the viscoelasticity moves into the bulk; The viscous drag significantly drops and the elastic wall 

Qi^ ' friction rises sharply. This transition is caused by an elastic force in the streamwise direction due 

to the competition of the convective and elastic effects. In the transition region, the convective and 
elastic effects balance. When the elastic effect dominates the convective effect, the second transition 
occurs but it is relatively moderate. The second one seems to be governed by so-called Weissenberg 
effect. These transitions are not sensitive to driving forces. By the scaling analysis, it is shown that 
the stress component is proportional to the Reynolds number on the boundary of the first transition 
in the Re-Wi space. This scaling coincides well with the numerical result. 

o 

^ : 1 Introduction 



Nonlinearities due to viscoelasticity in fluid flows lead to a variety of non-trivial flows such as rod- 
climbing, secondary flows, and tubeless siphon [I]. In polymeric fluids, even though the amount of the 

^^ ' polymer additives is very small, such a viscoelastic nonlinearity is effective. For decades, the under- 

standing of the formation mechanism of the non-trivial flows due to nonlinear viscoelasticity has been 
of academic and industrial interest. 

The nonlinearity originates from a viscoelastic stress augmented by the strain of flows and the 

/\^ • feedback of the elastic stress to flows. The target of our study is a flow pattern formation due to the 

nonlinear coupling in curvilinear channel shear flows. Flow pattern formations due to the nonlinear 
coupling are typically observed in rod climbing, i.e. Weissenberg effect [2j. The rod climbing is caused 
by an elastic stress working to the centripetal direction of curved streamlines (called "hoop stress") as 
follows. When a rotating rod is inserted into viscoelastic fluids, the rotating shaft forms shear flows in 
the fluids. Because the shear flows stretch and align polymer molecules along streamlines, an elastic 
stress works to the inward direction of streamlines by the tension of elongated polymer molecules. The 
fluid transferred to the inward direction finally climbs up along the rod. 

The earlier studies have shown that the hoop stress is also the primary effect inducing elastic insta- 
bilities in Taylor-Couette flows [31 13] . In this system, perturbed radial flows cause additional hoop stress 
and the hoop stress then further intensifies the radial flows. The elastic instabilities are followed by flow 
pattern formations including Taylor-Couette vortices, rotating standing waves, disordered oscillations, 
and time-dependent or time-invariant solitary vortices between the cylinders with rotational speed of an 



inner cylinder [5j |4j |6l [7] ■ Elastic instabilities in curvilinear shear flows are also observed in Taylor-Dean 
flow [5], cone-and-plate flow [S], and lid-driven cavity flow [TU] . 

These instabilities occur even at very low or even zero Reynolds number Re ^1. In this regime, 
elasticity-induced turbulence has also been observed in various experiments of shear flows including 
Taylor-Couette flow, von Karman swirling flow, and curvilinear channel flow [TTJ [T^l [13] ■ The character- 
istics of the elastic turbulence are an increased flow resistance and efflcient mixing at very low Reynolds 
number [14j . This mixing is very useful for industrial applications such as micro/nano devices [TF. 

Our goal of the investigations is the understanding of the detailed mechanism of flow pattern for- 
mation, elastic instability, and elastic turbulence in viscoelastic curvilinear channel flows as employed 
in the experiments of elastic turbulence [121 HH [H] • Although the previous studies concerning elastic 
instabilities in curved flows have basically noticed a radial force due to the hoop stress, in this study, 
we focus on the flow pattern transition caused by a streamwise elastic force, that is an elastic force in 
the streamwise direction in curved streamlines. The transitions are characterized by the disappearance 
of the separation vortex and the emergence of a jet flow with the sudden variation of wall frictions, and 
are caused by the increase of streamwise flows along curved streamlines due to the streamwise elastic 
force. Note that the flow pattern transition due to the streamwise elastic force is effective only when 
an inertial effect is also effective at moderate Reynolds number unlike the earlier elastic instabilities for 
Re <C 1. This indicates that the competition of the convective and elastic effects is essential for the flow 
pattern transition in the curvilinear channel shear flows. 

In order to investigate a detailed mechanism of flow pattern formation, direct numerical simulation 
is an effective approach on such complex flows. Nevertheless, earlier direct numerical simulations in 
the context of elastic instabilities have not treated such complex flows but treated more simple flows 
including two-dimensional Kolmogorov and four-mill flows driven by external forcing in doubly-periodic 
boundary conditions which have no boundary walls [71 1171 \TE\ . Therefore, we perform direct numerical 
simulations of periodic shear flows in two-dimensional wavy-walled channels as sketched in Fig. [I] as 
a model of curvilinear channel shear flows with high-accuracy spectral element method. The periodic 
wavy- walled geometry is canonical for flows with open and curved streamlines 19 and has actually also 
attracted much attention for industrial applications such as compact heat exchangers [201 121| . 

In this study, we investigate base laminar steady solutions for low to moderate Reynolds and Weis- 
senberg numbers. Note that Newtonian flows in this system is linearly unstable at moderate Reynolds 
number [5T] unlike two-dimensional Taylor-Couette and Taylor-Dean flows of viscoelastic fluids [^I23j . 
This linear instability may show that viscoelastic fluid flows are also linearly unstable in wavy-walled 
channels. We expect that the study of the base solutions leads to future elastic instability investigations. 

This paper is organized in the following way. Section[21provides governing equations and several non- 
dimensional parameters. A numerical method, boundary conditions and a flow geometry are described 
in Sec. [31 The validation of our computations is done and parameter values used in this study are 
summarized in Sec. [31 Section [3] shows the main results of our study. First, the steady solutions is 
classified into three groups according to the fiow patterns. Second, the variation of a viscous wall 
friction, a flow rate, and an elastic wall friction with Weissenberg number is described. Third, we 
describe the mechanism of the flow pattern transition and the variation of the above flow characteristics 
by utilizing the relationship among the norms of the terms in the Navier-Stokes equation and among 
local forces. Fourth, we perform scaling analyses in order to confirm the relationship of the norms. 
Fifth, the fiow pattern formation in Poiseuille-type flows driven by constant body force are also shown. 
Section [5] is devoted to discussion and concluding remarks. 

2 Basic equations 

The equations governing viscoelastic fluid flows consist of an equation of continuity, a momentum equa- 
tion and a constitutive equation describing the temporal evolution of a rheological state. We employ 
the incompressible Navier-Stokes equation for the momentum equation and the FENE-P model for the 



constitutive equation. The governing equations are written in the following dimensionless form: 
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where it is the velocity, p is the pressure, and t is the extra-stress tensor. C is the conformation tensor 
and L is the maximum polymer length. 

The equations include the following non-dimensional parameters: the Reynolds number, the Weis- 
senberg number, and the viscosity ratio parameter. The Reynolds number is the ratio of an inertial 
force to a viscous force defined in the following way: 



Re = pUoLo/ri, 



(6) 



where p is the fluid density, Lq is the typical length of flow geometry, and Uq is the typical flow speed. 
1] ^ rjg + rjp is the total viscosity defined as the sum of the solvent viscosity rjs and the polymer viscosity 
r]p. The Weissenberg number is the ratio of a polymer relaxation time and a typical flow time: 



Wi = XUo/Lo, 



(7) 



where A is the Zimm relaxation time. The viscosity ratio parameter is the ratio of the solvent viscosity 
and the total viscosity: 

/3 = Vs/V- (8) 

From the microscopic viewpoint, the FENE-P model is derived from a kinetic theory |i24j and the 
pre-averaging closure approximation introduced by Peterlin [3S] concerning elastic polymer dumbbells 
consisting of two beads and a spring connecting them. The conformation tensor C represents a measure 
of the second-order moment of the end-to-end distance vector of polymer dumbbells: 



Cij — {riTj), 



(9) 



where r is the relative difference of two beads positions. The function /(r) in Eq. ([5]) is so-called Peterlin 
function, which limits the length of polymers to be lower than the maximum polymer length L. When 
L ^ oo, the FENE-P model coincides with the Oldroyd-B model, which is derived from a kinetic theory 
on linear Hookean dumbbells. 

Because the FENE-P model reproduces the flow properties qualitatively similar to those observed 
in experiments, for example "shear thinning" |26j . the FENE-P model has been frequently employed in 
many studies on drag reduction [27l [28] , elastic instability [291 E] i and the other general viscoelastic 
fluid flows. We hence adopt the FENE-P model as a constitutive equation. 



2.1 Additional diffusive term 

In order to avoid numerical instabilities, we add an artificial diffusive term kV^C in the constitutive 
equation ([3]), where n is the artificial constitutive diffusivity coefficient. The artificial diffusive term 
added to the constitutive equation physically originates from the Brownian motion of the center-of- 
mass of elastic dumbbells across streamlines [30]. Sureshkumar [31] investigated the influence of the 
addition of the artificial diffusive term on flow stability in detail. They concluded that sufflciently small 
artificial diffusivity coefficient (k < 10^'^) does not affect the critical eigenmodes of the viscoelastic 
Orr-Sommcrfeld problem though the artificial diffusivity coefficient needed for numerical stabilization is 
much larger than the values physically obtained in the earlier study |30j . This implies that sufficiently 
small artificial diffusivity coefficient does not qualitatively change flow behavior of viscoelastic fluids 
from the viewpoint of numerical instabilities. This additional diffusivity hence has been widely used 
in the calculations of viscoelastic fluid flows governed with the Oldroyd-B, FENE-P, and other models 

[IH1I32]- 

3 Numerical method 

3.1 Spectral element methods 

We adopt spectral element method (hereinafter referred to as SEM) to numerically solve the governing 
equations in a complex geometry such as wavy- walled channels with high accuracy [331 IMl I3S]. In 
SEM, a computed domain is subdivided into several spectral elements. The values of variables and the 
derivatives of the variables are interpolated through the values on grid points with Lagrange polynomial 
interpolation functions in the elements [36] . Numerical errors between numerical solutions and accurate 
solutions are reduced in two ways in SEM. First, the increase of the polynomial degree of interpolation 
functions results in exponential convergence of numerical accuracy for sufficiently smooth solutions j37| . 
Second, the increase of the number of elements leads to an algebraic convergence similar to that of 
finite element calculations. SEM thereby enables the calculations with flexible boundary conditions 
and high numerical accuracy which are advantages of finite element methods and spectral methods, 
respectively. Because the local mesh refinement is required for viscoelastic calculations due to an elastic 
stress enhancement near walls, we provide more grid points near walls as shown in Fig. jljb). 

The Galerkin formulation of the primitive equations (IT])-(j3]) is used for spatial discretization in SEM. 
The Galerkin formulation in a two-dimensional domain 51 is written as follows: 

(V • u)qdVL = 0, (10) 
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where the tensor product is defined for two arbitrary tensors P and Q as: 

P■.Q = P^,Q^J, 1^1,3^2 (13) 
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Table 1: Grid point numbers and parameter values fixed in this study. A^e; is the number of elements 
and Ny is the order of interpolation polynomials of the velocity. 

and the notation (Vit)y = duj/dxj is adopted. The pressure p and the test function q belong to the 
functional space Q C L'^{Q). The velocity u and the test function w belong to the space W C ff^(n)^. 
The conformation tensor C and the test function S belong to S C 7?^ (17)^^^. Here L'^{i}) and H^{n) 
represent the space of square-integrable functions and the space of derivative functions whose first order 
derivatives are square-integrable on fl. The integral of these equations is performed with Lagrange 
interpolation functions and quadrature rules. Finally, nonlinear algebraic equations described with 
operator matrices and vectors of variables are derived. The linear terms are discretized by a second- 
order backward difference (BDF2) scheme and the nonlinear terms are discretized by a second-order 
extrapolation (EX2) scheme. We adopt a decoupled approach to solve time-splitting equations [55] . 

3.2 Boundary condition 

No-slip boundary conditions are imposed at the upper and bottom walls for the velocity. The boundary 
integral, which stems from the partial integral of the artificial diffusive term of the constitutive equation, 
assumes to be zero for numerical convenience. A periodic boundary condition is imposed on the velocity 
and the conformation tensor on infiow and outfiow boundaries. No boundary condition is imposed on 
the pressure because the pressure has no grid points on the boundaries. 

We mainly investigate Couette-type flows in which the upper wall is at rest and constant speed is 
given on the bottom wall. The wall speed of the bottom wall is normalized to unity along tangential 
directions t denoted in Fig. [TJa) at every point on the bottom wall. We also investigate Poiseuille-type 
flows in which both the upper and bottom walls are at rest. 

3.3 Parameters 

Here, we summarize the parameters used in this study: (1) fluid and rheological parameters are the 
Reynolds number Re, the Weissenberg number Wi, the viscosity ratio parameter /?, and the maximum 
polymer length L, (2) geometric parameters are the wall amplitude A and the periodic length L^ defined 
at Sec. 13.51 (3) the artificial diffusivity coefficient k. We focus on the dependence of fiow behavior on 
the Reynolds number, the Weissenberg number, and the curvature of streamlines controlled by the wall 
amplitude A. The other parameters /3, L, k, and L^ are fixed to the values shown in Tab. [TJ We set 
the maximum polymer length L to the sufficiently large value L^ = 1.0 x 10'' in order to capture the 
influence of a large elastic stress on fluid flows. The values of the fixed parameters /3, L, and k have been 
frequently used in the previous numerical studies of viscoelastic fluid flows [THl [Zl ISB 132 • Our several 
test runs show that these parameter values do not qualitatively influence the flow behavior obtained in 
this study. 

3.4 Validation of our computations 

Validation of our numerical code is done on plane-Couette flows. We perform the calculations at the 
low Weissenberg number Wi = 1.0 with zero artificial diffusivity coefficient k — and at the high 
Weissenberg number Wi = 20.0 with the finite artiflcial diffusivity coefficient k — 1.0 x lO^''. 

We compare our numerical solutions with exact steady solutions of the FENE-P model. The exact 
solutions of plane-Couette flows of the FENE-P model are derived as the roots of the following equations 
((T4)l and (fT5| obtained from the original equations under the condition that the derivatives of all variables 
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Table 2: Absolute value of the extra-stress component Tyy with several parameter values for validation 
of our numerical code. 

with respect to x are zero, p — const, u^ = \ — y^ and Uy — 0: 

2Wi2c3j^ + [L^K + 2)Cyy - L - 2 = 0, (14) 

C^y + WlCly = 0, a,^ ^L^ + (L^K + l)Cyy ^ 0. (15) 

The conformation tensor C of the exact solutions is uniform in the whole domain and the extra-stress 
component Tyy is zero everywhere. We calculate the extra-stress component Tyy for Validation of our 
computations. 

Table. [2] shows the absolute values of the extra-stress component Tyy with several parameter values. 
The difference between our results and the exact solutions is within 10^^ for the wide range of parameter 
values. Moreover, numerical accuracy is improved with mesh refinements. 

3.5 flow geometry 

The governing equations are solved in the domain sketched in Fig. [Ija). Two sinusoidal curves are 
boundary walls. The coordinates of the upper and bottom walls are defined as follows: 

(upper): y = 1 + Asin(27ra;/i2;), 
(bottom): y = As\n{2'Kx/ L^), 

where A is the amplitude of wall corrugations and L^ is the length of the periodic channel in the direction 
of X. The width of the channel is normalized to unity. We denote this channel "sinuous channel" as 
called by Cho f2T|. We call the regions denoted as "a" and "e" in Fig. [T](a) "outward region" because 
these regions are located in the outward position of curved flows parallel to the wall. We similarly call 
the regions denoted as " b" and " d" in Fig. [IJa) " inward region" . The regions denoted as " c" and " f ' in 
Fig.[IJa) are expressed "contraction region" in this paper. The symbols "t" and "n" denote a tangential 
and a normal directions to the wall, respectively. 

We adopt the sinuous channel in this study for the following reasons. First, the channel has mixed- 
kinematics with both shear and extensional flows. The elastic stress due to polymer stretching in 
these flows strongly influences flow behavior. Second, we can easily control the magnitude of shear 
and extension by changing the wall amplitude A and the periodic length L^ of the channel. Third, 
streamlines are highly curved in the channel. The curvature of streamlines locally varies and the sign 
of the curvature is even reversed along the streamlines unlike two-dimensional Taylor-Couette flows in 
which the sign of the curvature is not changed along streamlines. The spatial variation of the curvature 
leads to the change of the magnitude and the direction of an elastic hoop stress referred at Sec. [T] The 
spatial variation of the hoop stress plays a crucial role in the formation of flow patterns as described 
below. 

Despite the importance of the sinuous channel, most of the earlier works of viscoclastic fluid flows 
in wavy-walled channels have considered corrugated (contraction/expansion) channel flows in which the 
phases of sinusoidal curves of an upper and a bottom walls are shifted in a half wavelength [39l [191 SOI HI] ■ 
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Figure 1: (Color online) (a) Flow domain with symbols and expressions. The regions "a" and "e" are 
the outward regions, "b" and "d" are the inward regions, and "c" and "f are the contraction regions. 
The symbols "t" and "n" denote a tangential and a normal directions to the boundary wall, respectively, 
(b) Computational mesh used in this study. 

The studies of the sinuous channel flows are limited to those of Newtonian fluid flows [501 [51] ■ However, 
the corrugated channel is not appropriate for our study focusing on the curvilinear channel flows because 
the streamlines of a bulk flow in the corrugated channel are more weakly curved than those in the sinuous 
channel. 



4 Results 



4.1 Flow pattern transition and classification of steady solutions 

We perform direct numerical simulations for the wide range of Weissenberg number < Wi < 20 and 
Reynolds number < Re < 1000. We obtain steady solutions in these regions. 

The steady solutions are classified into three groups. The three groups correspond to three regions in 
the Wi-Re parameter space as shown in Fig.[5| We call the three regions " convective" , "transition", and 
"elastic" regions. The two Weissenberg numbers separating adjacent regions is called first and second 
critical Weissenberg numbers, Wici and Wic2, defined at Sec. 14.31 Wic2 is extrapolated into a large- Re 
region using a fitting function because our numerical calculations suffer numerical instabilities at large 
Wi and Re shown in the region with diagonal lines in Fig. [5] 

The difference of the steady solutions in the three groups typically appears in the fiow patterns of 
the steady solutions. Figures [3| and [4| show streamlines and flow-speed profiles of the steady solutions 
in the convective region Wi = 1.0, the transition region Wi = 5.16, and the elastic region Wi = 8.0 at 
the moderate Reynolds numbers Re = 200 and 400. Figures [3K a) and[3][a) show that, in the convective 
region, the fiow field consists of a bulk shear flow and a vortex in the outward region (a) (see the 
definition of the outward region (a) in Fig. [IJa)). We call the vortex "separation vortex" because the 
vortex is detached from the bulk flow in the channel. The separation vortex has a symmetrical shape 
at low Wi, while the vortex becomes smaller and the symmetry of the vortex becomes gradually broken 
with Wi. The separation vortex is finally vanished or becomes very small near the upper wall at Wid 
as shown in Fig. [3Kb). As discussed in detail, the depression of the separation vortex is just caused by 
the increase of fiow due to the streamwise elastic force in the outward region (a). The separation vortex 
reappears and subsequently grows with Wi as shown in Fig.[3Kc). 

The other difference of fiow patterns in the three regions is "jet flow" observed in the vertical section 
of flow velocity. Figures [Sj [6l and [7| show the vertical sections of the velocity component Ux at x = 7r/4 
with Re = 200, 400, and 800. We define jet fiows the region in which the velocity is larger than in 
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Figure 2: (Color online) Phase diagram of the base steady solutions in the plane of Wi and Re. The left 
and right curves denote the first and second critical Weissenberg numbers, Wici and Wic2, defined at 
Sec. 14.31 respectively. The Reynolds number of the edge point on the left curve Re = 70 is determined 
through the scaling analysis at Sec. 14.41 The curve of Wic2 is extrapolated to a large-Re region using 
a fitting function. The value of Wic2 at Re = is the Weissenberg number at which the separation 
vortex emerges. The region with diagonal lines represents that in which our calculations suffer numerical 
instabilities. 



the surrounding area. Jet flows whose velocity is larger than wall speed (unity) are observed near the 
bottom wall 0.5 < y < 0.7 in the convective region as shown in the curves [a]-[c] in Figs. [SJa), El^a), 
and[7Ja). The jet flows near the bottom wall suddenly disappear at Wid as shown in the curve [d] in 
Figs.[Sl[a),[nia), and[71[a). Simultaneously, jet flows appear near the upper wall as shown in the same 
velocity profiles. The jet flows near the upper wall are broader than those near the bottom wall in the 
convective region. The jet flows gradually reach to the bottom wall with Wi during the transition region 
as shown in the curves [d] and [e] in Figs. [SJa), EJa), andjTJa). In the elastic region, the jet flows are 
repeatedly present near the bottom wall though the jet flows are slower and broader than those in the 
convective region as shown in the curve [f] in Fig. Efa), [St^a), and[71[a). The emergence of the jet flows 
does not depend on Re except the jet flows at larger Re is broader in the elastic region. As described 
below in detail, the jet flows are also caused by an streamwise elastic force due to polymer stretching 
in viscoelastic fluid flows. Note that negative-velocity regions appear near the upper wall as shown in 
the velocity profiles magnified near the upper wall in Figs.[5jb),|6Ub), and [71(b). These negative- velocity 
regions exist in the separation vortices in the outward region (a) as shown in Figs.[3]and|4l 

The jet flows are also observed in the flow-speed profiles as shown in the right figures in Figs. |3]and 
m In the convective region, the bright regions in the inward region (b) correspond to the jet flows. In 
the transition region, the jet flows correspond to the regions where the flow speed is larger than the 
surrounding area, from the outward region (a) to the contraction region (c) near the upper wall. In the 
elastic region, the jet flows emerge from the contraction region (c) to the contraction region (f) via the 
inward region (b) near the bottom wall. 

We call the change of flow patterns including the separation vortex and the jet flow in the three 
regions " flow pattern transition" . 
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Figure 3: (Color online) Streamlines (left) and flow-speed profiles (right) of the steady solutions in (a) 
the convective region Wi — 1.0, (b) the transition region Wi = 5.16, and (c) the elastic region Wi — 8.0 
at Re — 200. Streamlines are drawn using passively-advected particles. The intervals of streamlines in 
bulk flows are determined to be inversely proportional to the flow speed at the position. The separation 
vortices are drawn independently of the bulk flows. 
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Figure 4: (Color online) Streamlines (left) and flow-speed profiles (right) of the steady solutions in (a) 
the convective region Wi = 1.0, (b) the transition region Wi — 5.6, and (c) the elastic region Wi — 8.0 
at Re — 400. The streamlines are drawn with the same method in Fig. [31 
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[a] Wi=1.0 

[b] Wi=4.0 
c]Wi=5.14 

[d]Wi=5.16 
Wi=6.0 
[f] Wi=8.0 
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[a] Wi=1.0 

[b] Wi=4.0 
[c]Wi=5.14 
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Figure 5: (Color online) (a) Tangential velocities ut in the vertical section at cc = 7r/4 in the convective 
region [a] Wi = 1.0 [b] Wi = 4.0, [c] Wi = 5.1, in the transition region [d] Wi = 5.16, [e] Wi = 5.5, 
and in the elastic region [f] Wi — 8.0 at Re — 200. (b) Tangential velocities Ut in the vertical section 
at X = 7r/4 magnified in 1.3 < y < 1.5. (c) Tangential velocities ut in the elastic region [a],[b] Wi = 8.0 
and [c],[d] Wi = 10.0 in the vertical section at [a],[c] x = 7r/4 and [b],[d] x = 7r/2. 
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Figure 6: (Color online) (a) Tangential velocities ut in the vertical section at a; = 7r/4 in the convective 
region [a] Wi = 1.0 [b] Wi = 4.0, [c] Wi — 5.5, in the transition region [d] Wi — 5.6, [f] Wi = 6.0, and 
in the elastic region [e] Wi = 10.0 at Re = 400. (b) Tangential velocities ut in the vertical section at 
X = 7r/4 magnified in 1.15 < y < 1.5. 
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Figure 7: (Color online) (a) Tangential velocities Ut in the vertical section at a; = tt/4 in the convective 
region [a] Wi = 1.0 [b] Wi = 4.0, [c] Wi = 6.04, in the transition region [d] Wi = 6.06, and [e] Wi = 7.0 
at Re = 800. (b) Tangential velocities ut in the vertical section at x = 7r/4 magnified in 1.2 < y < 1.5. 
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Figure 8: (Color online) Average viscous wall friction r^^u/ as a function of Wi at Re = 200 and 400. 
Two vertical lines represent the first and second critical Weissenberg numbers at Re = 200. 

4.2 Characteristics of flows 

In the previous subsection, we described that the steady solutions are classified into the three groups. 
The solutions belonging to different groups have remarkably different flow patterns. How the change of 
flow patterns influences characteristic quantities of flows such as a wall friction and a flow rate? In this 
subsection, we show the dependence of a viscous wall friction, an elastic wall friction, and a flow rate 
on Wi at Re = 200 and 400. 

4.2.1 Laminar drag reduction 

The flow pattern transition described in the previous subsection accompanies the abrupt change of a 
viscous wall friction. The viscous wall friction is very important for engineering applications to the 
control of the loss of flow rate. 

We define the average viscous wall friction (viscous drag) acting on the walls of the channel as follows: 



1 



dut 



dn 



dr, 



(16) 



where F represents the total of the bottom boundary Fi and the upper boundary F2. Sr is the sum 
of the lengths of the upper and bottom boundary walls, t and n denote a tangential and a normal 
directions as shown in Fig. [TJa). dut/dn represent the derivative of the tangential velocity component 
Ut with respect to the wall-normal direction n. 

Figure |S] shows the average viscous wall friction as a function of Wi at Re = 200 and 400. The 
viscous drag is abruptly reduced at Wici, while the viscous drag gradually grows in the convective and 
elastic regions. The abrupt reduction of the viscous drag simultaneously occurs with the vanishment of 
the separation vortex. 

We call the drag reduction in laminar steady fiows " laminar drag reduction" . The expression has 
been also used in the studies of viscoelastic fluid flows in curved pipes with circular cross sections [12] . 
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Figure 9: (Color online) Flow rate FR as a function of Wi at Re — 200 and 400. Two vertical lines 
represent the first and second critical Weissenberg numbers at Re = 200. 

4.2.2 Flow rate 

The flow rate FR in the channel is defined in the following way: 



FR 



v^{7r/2, y)dy. 



(17) 



Figure [3] shows the flow rate FR as a function of Wi at Re — 200 and 400. In the convective region, 
the flow rate gradually increases, while the flow rate is temporarily decreased at Wici. The flow rate 
greatly grows in the transition region. In the elastic region, the growth of the flow rate is relaxed. 



4.2.3 Growth of flow resistance 

In this subsection, we evaluate the following average elastic wall friction: 



"^ewf 



i/r"*"'^' 



(18) 



where rt„ is a component of the extra-stress tensor in the coordinate which consists of two base vectors. 
One is a tangential unit vector (denoted as t) and the other is a normal unit vector (denoted as n) as 
shown in Fig.[IJa). This elastic wall friction represents a flow resistance resulting from the elastic stress. 
The sum of the viscous and elastic wall frictions is a total flow resistance acting on the walls. 

Figure [To] shows the elastic wall friction as a function of Wi at Re = 200 and 400. At Wici, the 
elastic wall friction suddenly grows. The sudden growth of the flow resistance also simultaneously occurs 
with the vanishment of the separation vortex described at Sec. 14.11 

4.3 Contribution of elastic force onto the flow pattern transition 

In this subsection, we first introduce the norm of the absolute value of each term in the Navier-Stokes 
equation. We next elucidate the mechanism of the flow pattern transition by utilizing the relation among 
the norms. 

The norm of the absolute value of each term is defined in the following way: 



|Re(tt- V)m|| 



\Re{u-V)u\dn. 



(19) 
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Figure 10: (Color online) Average elastic wall friction Te^,/ as a function of Wi at Re = 200 and 400. 
Two vertical lines represent the first and second critical Weissenberg numbers at Re = 200. 

Figure [TT] shows the absolute- value norms of the convective term ||Re(M- V)m||, the pressure gradient 
term ||Vp||, the elastic term |1V • t||, and the viscous term ||/3V^m|1 in the Navier-Stokes equation as a 
function of Wi at Re = 200. 

We define the first critical Weissenberg number Wici as the Weissenberg number at which the norm 
of the elastic term intersects with that of the pressure gradient term. The second critical Weissenberg 
number Wic2 is similarly defined as the Weissenberg number at which the norm of the pressure gradient 
term crosses with that of the convective term. 

We calculate the maximum speed of flows in the separation vortex as an intensity of the separation 
vortex. The maximum speed is defined in the following way. We consider the region where Ua; < in 
the channel. This region exists in the separation vortex as shown in Figs.jHJa) andlS^b). We calculate 
the maximum value of the flow speed \u\ in this region. When the separation vortex is vanished, the 
region where u^; < also disappears. Then, the maximum speed is set to zero. In Fig. 1121 we show 
the maximum speed as a function of Wi at Re = 100 and 200. This figure shows that the first critical 
Weissenberg number Wici coincides with the Weissenberg number at which the separation is suddenly 
depressed. 

The second critical Weissenberg number provides an estimate of the Weissenberg number at which 
the jet flows reach to the bottom wall in the transition region, though the estimate is less exact than 
the accordance of Wici to the Weissenberg number at which the separation vortex is depressed. 

We subsequently describe the mechanism of the flow pattern transition. Here, we utilize the profiles 
of the pressure and the conformation tensor as well as the norms of the terms for the description. Figures 
[T51 and [m display the profiles of the pressure and the trace of the conformation tensor Tr(C) in the 
convective, transition, and elastic regions at Re = 200. The trace of the conformation tensor represents 
the magnitude of polymer stretching. Note that the pressure is fixed to zero at the given point. Negative 
pressure at certain point then represents the pressure smaller than at the reference point. In these 
figures, we can see the local relation among elastic force, pressure gradient, and centrifugal force defined 
as Re(U( /i?) where R is the curvature radius of streamlines. Here, we notice that the streamwise elastic 
force plays a key role in the flow pattern transition as well as the hoop stress. Figure [TS] schematically 
summarizes the transition of the local relation among the forces in the three regions. At the subsequent 
subsections, we describe the flow pattern transition in the three regions by utilizing these quantities. 
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4.3.1 Convective region 

In the convective region, the separation vortex emerges in the outward region (a) (see Fig.lHJa)). This 
is a characteristic phenomenon of Newtonian fluid flows in wavy- walled channels |20j . 

The generation mechanism of the separation vortex is as follows. The norms of the convective and 
pressure-gradient terms almost accord in the convective region as shown in Fig. llll while the elastic and 
viscous terms are much less than the convective and pressure- gradient terms. This accordance indicates 
that the pressure gradient and the centrifugal force cancel in the whole domain in steady flows. The 
centrifugal force is very large particularly in the inward region (b) where the streamwise velocity and 
the curvature are large. The large centrifugal force leads to the large centripetal pressure gradient in 
the inward region (b) . Then, the pressure becomes smaller in the inward region (b) than in the outward 
region (a) (see Fig. [T5Fa)). Because the pressure gradient inhibits the flows from the contraction region 
(f) to the outward region (a), the bulk flow then passes near the inward region (b). The separation 
vortex detached from the bulk flow finally emerges in the outward region (a) (see Fig. [3][a)). 

In the convective region, the characteristic phenomenon of viscoelastic fluid flows is observed as 
a jet flow near the bottom wall as shown in Figs. [5ja), EJa), and[71[a). The generation mechanism 
of the jet flow is described as follows. The bulk shear flow through the inward region (b) strongly 
stretches polymers near the bottom wall. The polymer stretching leads to the increase of the elastic 
stress and to that of the streamwise elastic force in this region (see Fig. HU^a)). Because the augmented 
streamwise elastic force increases a streamwise flow, the jet flow finally emerges in the inward region (b) 
(see 0.5 < y < 0.7 in the curves [a]-[c] in Fig. [SJa)). 

The jet flow causes the viscous drag enhancement in the convective region as described at Sec. 14.2.11 
in the following way. Due to the jet flow, the wall-normal velocity gradient is inverted near the bottom 
wall. The elastic force near the bottom wall gradually increases with Wi (see Fig. fTU)) . Then, the velocity 
of the jet flow is also augmented due to the intensified elastic force as shown in the curves [a] and [b] 
in Fig. Elja). Because the increase of the velocity leads to that of the velocity gradient, the viscous wall 
friction also grows on the bottom wall as shown in Fig. |8l The increase of the flow velocity due to the 
streamwise elastic force near the bottom wall also directly leads to the growth of the flow rate as shown 
in Fig. M 

With further increasing Wi, the flow behavior described above is moderately modified. The repre- 
sentative modification is the growth of the elastic stress in the contraction regions (c) and (f ) (see Fig. 
I14r b)). In the contraction regions (c) and (f), the channel is narrower than in the highly-curved regions 
including the outward regions (a), (e), the inward regions (b), and (d) (see Fig. [ija)). The flow velocity 
is basically large in the narrow regions (c) and (f) because the flow rate, which is the integral of the 
flow velocity in the wall-normal section of the channel, is constant in the every wall-normal section. 
The difference between the velocities in the contraction region and in the highly-curved region stretches 
polymers between the two regions. The augmented elastic stress forms the lumps of large-stress regions 
in the contraction regions (c) and (f) (see Fig. ll4r b)). The lump in the contraction region (c) gradually 
extends to the outward regions (a) and (e) with Wi and the lump in the contraction region (f ) gradually 
reaches to the inward region (b) (see Fig. [TWb)). Because the lump of the elastic stress induces flows 
from the outward region (a) to the contraction region (c) due to the tangential gradient of the elastic 
stress (see the definition of the tangential direction in Fig. [Ija)), the separation vortex in the outward 
region (a) gradually becomes small with Wi (see the curves [b] and [c] in Fig.[5jb)). Since the large-stress 
region extended from the contraction region (f) to the inward region (b) increases flows above the jet 
flow in the inward region (b), the jet flow is relatively depressed (see 0.5 < y < 0.7 in the curves [b] and 
[c] in Fig. EJa)). The depression of the jet flow directly leads to the decrease of the viscous wall friction 
on the bottom wall at Wi > 4.0 as shown in Fig. [51 

The description at this subsection is summarized as follows (see also Fig. [TSja)). The separation 
vortex in the outward region (a) and the jet flow in the inward region (b) are the characteristic flow 
patterns of the steady solutions in the convective region. The separation vortex is caused by the local 
forces which also works in the Newtonian fluid flows, whereas the jet flow stems from the elastic force 
characteristic of viscoelastic fluid flows. 
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4.3.2 Transition region 

As described at the previous subsection, ttie flow from tiie outward region (a) to the contraction region 
(c) gradually increases with Wi due to an streamwise elastic force. The increased flow further augments 
the elastic stress in the contraction region (c) because the increased velocity and velocity gradient further 
stretch polymers. The augmented elastic stress repeatedly increases the flow from the outward region 
(a) to the contraction region (c) . This argument indicates that the flow and the streamwise elastic force 
is self-enhanced in this region. 

When Wi attains the first critical Weissenberg number Wici, this self-enhancing procedure actually 
sets in. Then, the elastic effect becomes suddenly comparable with the convective effect, and the flow 
pattern observed in the convective region drastically changes. 

The sudden formation of the jet flow in the outward region (a) and the sudden vanishment of the 
jet flow in the inward region (b) are the prominent change of the flow pattern in the transition region 
as described at Sec. 14.11 (see 0.5 < y < 0.7 and 1.1 < y < 1.5 in the curves [c] and [d] in Fig. [DJa)). 
The generation of the jet flow in the outward region (a) results from the increase of the velocity in the 
outward region (a) due to the streamwise elastic force described above. Then, the jet flow in the inward 
region (b) is relatively depressed due to the increase of flows to the outward region (a). The vanishment 
of the jet flow in the inward region (b) leads to the abrupt reduction of the viscous drag on the bottom 
wall as shown in Fig. |S1 Due to the emergence of the jet flow in the outward region (a), the separation 
vortex is vanished or becomes very small as shown in Fig.[21[b) and y > 1.4 in the curve [d] in Fig.jSJa). 

The lump of the large-stress region in the contraction region (c) also suddenly grows beyond Wici. 
Because the lump of the large-stress region attains to the walls (see Fig. [Ml c)). the elastic wall friction 
suddenly grows at Wici as shown in Fig. [TUl Note that the lump of the elastic stress is inclined to the 
channel walls. This is because the jet flow in the outward region (a) enhances the elastic stress near the 
upper wall (see Figs.fTJTc)). 

The flow rate for Wi > 4.0 varies with the decrease of the flow velocity in the inward region (b) 
and the increase of the flow velocity in the outward region (a) . The flow pattern transition leads to the 
reduction of the flow rate as shown in Fig. [HI However, the flow rate remarkably increases above Wici 
because the streamwise elastic force further increases the flow to the outward region (a) with Wi in the 
transition region. The growth of the flow rate with Wi also appears in the amplitude of the jet flow as 
shown in the curves [d]-[f] in Fig. [5ja). 

In the transition region, the jet flow in the outward region (a) gradually moves from the upper 
wall to the bottom wall (see the curves [d]-[f] in Fig. [5ja)). The mechanism of the moving jet flow is 
described as follows. The hoop stress induces the elastic force working to the centripetal direction in 
curved streamlines. Due to the centripetal elastic force in the outward region (a), flows normal to the 
streamwise direction are induced. Because of the flow normal to the streamwise direction, the lump 
of the large-stress region in the contraction region (c) falls down to the bottom wall (see Fig. EKd)). 
The descended lump of the large-stress region repeatedly induces the jet flow due to the streamwise 
elastic force in the lower position of the channel (see Fig. [T^d) and the curve [e] in Fig. EJa)). Because 
the jet flow augments the elastic stress between the outward region (b) and the inward region (a), the 
large-stress belt-like region is flnally formed from the region between the outward region (b) and the 
inward region (a) to the inward region (d) (see Fig.jTU^d)). The hoop stress strongly works in this region 
due to the large elastic stress and the large flow curvature. Then, the flow normal to the streamwise 
direction is further induced. Thus, the jet flow gradually reaches to the bottom wall. With the variation 
of the position of the jet flow, the separation vortex repeatedly emerges and gradually grows (see the 
curves [d]-[f] in Fig. [5Kb)). 

The influence of the elastic stress on the flow is confirmed through the comparison of the elastic force 
with the centrifugal force and the pressure gradient in the following way (see also Fig. fTSl b')'). First, the 
hoop stress in the outward region (a) described above is canceled by the centrifugal pressure gradient 
in the steady fiows because the centrifugal force is very small due to the small velocity near the upper 
wall. The pressure distribution shown in Fig. I13f b) indeed indicates the centrifugal pressure gradient 
in the outward region (a). Second, in the inward region (b), the sum of the centripetal elastic force and 
pressure gradient cancels the centrifugal force in the steady fiows because the convective term, the origin 
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of the centrifugal force, is yet larger than the elastic and pressure gradient terms at relatively small Wi 
in the transition region as shown in Fig. [TT] Since the centrifugal force does not almost depends on Wi 
(see Fig. [TT]), the growth of the elastic force directly leads to the decrease of the pressure gradient in 
the inward region (b) (see Figs. I13f a) and (b)). The reduction of the pressure gradient in the inward 
region (b) is also consistent with the decrease of the flow to the inward region (b) . 

The description at this subsection is summarized as follows (see also Fig. [TSl b)). The sudden 
vanishment of the jet flow in the inward region (b) and the abrupt formation of the jet flow in the 
outward region (a) due to the streamwise elastic force simultaneously occurs at the beginning of the 
transition region. The jet flow gradually reaches from the upper wall to the bottom wall due to the 
elastic hoop stress in the transition region. 

4.3.3 Elastic region 

The jet flow approaching to the bottom wall stops approximately above the second critical Weissenberg 
number. In the elastic region, the jet flow near the bottom wall and the separation vortex in the outward 
region (a) characterizes the flow patterns as in the convective region (see Figs.lHKa) and (c)), 

Despite the similarity of the flow patterns between the convective and elastic regions, the generation 
mechanism of the separation vortex and the jet flow is remarkably different in the two regions. In the 
elastic region, the large-stress belt-like region described at the previous subsection covers wide regions 
including the inward regions (b), (d) , the contraction regions (c), and (f) in the channel due to the large 
elastic hoop stress at large Wi (see Fig. [T3Ke)). Because the convective term does not vary with Wi as 
described above, the pressure gradient term grows with the elastic stress in the steady flows in order to 
cancel the elastic term. The large-pressure region thus also spreads into the wide regions of the domain 
(see Fig. [TBTc')'). Then, the flow from the outward region (a) to the contraction region (c) is severely 
inhibited by the large pressure gradient from the contraction region (c) to the outward region (a). Since 
the bulk flow cannot pass through the outward region (a) , the separation vortex is stably maintained in 
the outward region (a) . Note that the separation vortex appears in the low-pressure region in the elastic 
region unlike the separation vortex observed in the high-pressure region in the convective region (see 
Figs.[3l[a),(c), Figs. flSl a'). and (c)). The vortex emergence in the low-pressure region is a characteristic 
phenomenon of viscoelastic fluid flows. 

The difference of flows between the convective and elastic regions also arises in the viscous wall 
friction and the flow rate (see Figs. |S]and|ni). In the elastic region, the flow is basically developed in 
the inward regions due to the large elastic hoop stress. The larger streamwise elastic force at larger Wi 
causes the larger velocity of the jet flow and the larger velocity gradient in the inward region (b) (see 
Fig. EJc)). The variation of the velocity gradient directly leads to the viscous drag enhancement with 
Wi in the elastic region as shown in Fig. [51 Because the velocity gradient in the inward region (b) in the 
elastic region is much smaller than that in the convective region, the magnitude of the viscous drag in 
the elastic region is also much smaller than that in the convective region (see the curves [a]-[c] and [f] in 
Fig.[5][a)). The flow rates in the convective and elastic regions are also different because the magnitudes 
of the elastic stress in the two regions are remarkably different (see Figs.[T4Ua) and (e)). However, the 
growth rate of the flow rate is relaxed with Wi in the elastic region. This is because the elastic stress is 
also relaxed for sufficiently large Wi in the elastic region as shown in Fig. [TUj 

To summarize, the separation vortex in the low-pressure region and the jet flow in the inward regions 
(b) emerge due to the large elastic hoop stress and the streamwise elastic force in the elastic region (see 
also Fig. [TStc)). 

We have described the mechanism of the flow pattern transition, the variation of the wall frictions, 
and that of the flow rate in the three regions. These arguments show that the streamwise elastic force 
produces jet flows which characterize the flow pattern formation in wavy-walled channel shear flows. 

4.4 Scaling analysis on the elastic wall friction 

How is the relationship between the flow pattern transition and the change of the local force balance 
theoretically conflrmed? We perform several scaling analyses on the basis of the results described above 
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Figure 16: (Color online) First critical Weissenberg number Wici as a function of Re. The curve of Wici 
is fitted with a power function Fi(Re) == yli(Re - Rcc)^" + A3 where Rcc = 70, A^ ^ 0.06, Ai ^ 0.5 
and A3 ^ 4.5. 

for theoretical evaluation because it is severely difficult to derive analytical solutions in our flow system. 

In Fig. [121 we show the first critical Weissenberg number Wici as a function of Re. The curve of Wici 
is fitted by using a power function Fi(Re) = Ai(Re — Rcc)^^ + A3, where Rec is fixed to 70. The values 
of the fitting parameters are Ai ~ 0.06, A2 ^ 0.5, and A3 ~ 4.5. The relationship Wici ex \J^& — Rcc 
is nearly satisfied. Rec = 70 is a reasonable value because the sudden growth of the elastic wall friction 
does not occur when Re < Rcc . Moreover the size of the separation vortex in the convective region is 
abruptly decreased when Re < RCc. Although we draw a line from Re = 70 to Re = in a low Wi and 
low Re region in Fig. [2l the Reynolds number below which the separation vortex completely disappears 
in the Newtonian limit Wi — )• is maybe larger than zero. 

The norms of the convective term ||Re(M • V)tt|j and the elastic term |V • t\ have similar values at 
Wici as shown in Fig. [TT] When the non-dimensional velocity and the derivatives are scaled to 0(1), 
the following relationship is derived: 

T ex Re. (20) 

Figure [IT] shows the elastic wall friction rem/ at Wici as a function of Re. We superpose the linear line 
F2 = Si(Re-Rec) + B2 where Rcc = 70, B^ = 4.1 x 10"^, B2 = 3.6 on Fig.[T71 The linear relationship 
(PO)) between r and Re is indeed approximately satisfied. The two relationships r ex Re — Rcc and 
Wici ex Re — Rcc lead to the relationship of the elastic stress and Wici as t ex Wici^. 

4.5 purely elastic limit 

In this subsection, we show the results of purely elastic calculations at Re = 0. Recent years, purely 
elastic flows have been intensively studied in the context of elastic instability and elastic turbulence 
[T51 [TUl 131 HH] . The elastic force mainly drives flows in the purely elastic limit because the convective 
term is vanished. 

Figure [T51 shows streamlines of the steady solutions at the low Weissenberg number Wi = 1.0 and 
the moderate Weissenberg number Wi = 6.0 with Re = 0. The separation vortex is not observed at 
low Wi. The disappearance of the separation vortex results from the vanishment of the convective 
term. When the convective term is vanished, the centrifugal force does not work in the whole domain. 
Because the centrifugal force is canceled by the pressure gradient for low Wi as described at Sec. 14.3.11 
the pressure gradient is also remarkably weakened in the whole domain. The depression of the pressure 
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Figure 17: (Color online) Elastic wall friction Tewf at Wid as a function of Re. The linear line F2 = 
Bi{Re - Rcc) + B2 where Rcc = 70, Bi = 4.1 x 10"^ B2 = 3.6 is superposed. 

gradient which inhibits the flows from the contraction region (f) to the outward region (a) leads to the 
vanishment of the separation vortex in the outward region (a) at low Wi. 

The Weissenberg number between the transition and elastic regions at Re = in Fig. [3] nearly 
corresponds to the Weissenberg number at which the separation vortex appears. The flows at Re = 
are similar to those at moderate Re for large Wi as shown in Figs.jSjc) and lTST b). This is because the 
elastic force and the pressure gradient are dominant for large Wi at every Re. 

Figure [TOl shows the elastic wall friction as a function of Wi. The sudden growth of the elastic wall 
friction is not observed because the flow pattern transition between the convective and transition regions 
does not occur. However, the growth rates of the elastic wall frictions at the low- and high- Wi regions 
are much different. The Weissenberg number at which the growth rate changes nearly corresponds to 
the Weissenberg number at which the separation vortex emerges. 

4.6 Poiseuille-type flows 

In this subsection, we show the steady solutions of Poiseuille-type flows driven by a body force in order 
to elucidate the dependence of the flow pattern transition on driving forces. 

We add a constant body force (constant pressure gradient) only in the rectangular regions connecting 
the inflow and outflow boundaries. No body force is added in the two outward regions (a) and (e) (see 
Figs. [TJa) and 1^0)) . The wall amplitude A needs to be less than the half width of the channel A = 0.5 
in this method. We adopt A = 0.3 in the Poiseuille-type flow cases instead oi A = 0.5 used in the 
Couette-type flow cases. The magnitude of the body force is determined under the condition that the 
flow rate FR is fixed at every Wi. 

Figure [21] shows streamlines of the steady solutions of the Poiseuille-type flows at Re = 200. The two 
separation vortices appear in the outward regions (a) and (e) at the low Weissenberg number Wi = 1.0. 
These separation vortices disappear at the moderate Weissenberg number Wi — 6.0, but reappear at 
the high Weissenberg number Wi = 30.0. This shows that the flow pattern transition is not influenced 
by the difference of the driving forces between the Couette-type and Poiseuille-type flows. 

Next, we confirm the relationship between the fiow pattern transition and the variation of the elastic 
wall friction in the Poiseuille-type flows. In Fig. [22j we show the elastic wall friction as a function of 
Wi. When the separation vortices disappear, the growth rate of the elastic wall friction suddenly varies 
at Wi = 6.0. This result is also qualitatively consistent with that of the Couette-type flows. 
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Figure 18: (Color online) Streamlines of the steady solutions at the small Weissenberg number Wi = 1.0 
and the large Weissenberg number Wi = 6.0 with Re = 0. 
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Figure 19: (Color online) Elastic wall friction Tewf as a function of Wi at Re = 0. 



25 




Figure 20: (Color online) Regions in which a constant body force is imposed in Poiseuille-type flow 
cases. 

Although the Poiseuille-type and Couette-type flows are qualitatively similar, the Weissenberg num- 
bers at which the separation vortices disappear or emerge are much different between the two flows. The 
difference may result from that of the relationship between the centrifugal force and the elastic force. In 
Poiseuille-type flows, the centrifugal force is larger in a bulk region than near walls and the elastic force 
is larger near walls than in the bulk region. In contrast, the centrifugal force is larger near the bottom 
wall than in a bulk region and the elastic force is large in the inward regions (b) , (d) , in the contraction 
regions (c), and (f) in Couette-type flows. 

5 Discussion and Concluding remarks 

In this study, we seek base steady solutions of the curvilinear shear flows in the wavy-walled channel. 
However, the transient behavior of the relaxation process to the steady solutions is also important when 
we consider the relation of the base solutions to elastic instability. In this section, we first discuss the 
transient behavior using relaxation times against small perturbations on the base steady solutions and 
next summarize our study. 

In Fig. 1231 we show the relaxation processes of the velocity component u^ against small random 
perturbations on the conformation tensor of the base steady solutions in the convective region Wi = 2.0, 
the transition region Wi = 5.2, and the elastic region Wi = 8.0 at Re = 200. The same perturbation is 
added to the three steady solutions. The time is normalized with Wi which represents the relaxation time 
of the conformation tensor (see Eqs. ([3]) and ([7])). This figure shows that the perturbation more rapidly 
relaxes in the convective and elastic regions than in the transition region. Note that the perturbation 
amplitude of the velocity is the largest in the transition region despite the same amplitude of the stress 
perturbations. The large velocity perturbations shows that the feedback from the elastic stress to the 
velocity is the largest in the transition region. We expect that the pressure gradient plays a key role 
in the difference of the relaxation processes between the three regions. In the convective and elastic 
regions, the pressure gradient nearly accords with the centrifugal force and the elastic force, respectively. 
In contrast, the pressure gradient in the transition region relatively decreases compared with the other 
regions (see Fig. fTTj) . The pressure reduction and the long relaxation time in the transition region may 
indicate that the pressure gradient works to depress the perturbation of the extra-stress in our flow 
system. Note that the most oscillatory behavior with the smallest period is observed in the elastic 



26 




Figure 21: (Color online) Streamlines of the Poiseuille-type flows at the low Weissenberg number Wi = 
1.0, the moderate Weissenberg number Wi — 6.0, and the high Weissenberg number Wi = 30.0 at 
Re = 200. 
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Figure 22: (Color online) Elastic wall friction r^wf in the Poiseuille-type flows as a function of Wi. 

region as shown in the inset of Fig. EHT c). The oscillatory nature in the high-Wi region is also observed 
in the earlier two-dimensional calculations of viscoelastic extensional flows |18j . 

In order to sustain the disturbance against the relaxation due to the pressure gradient, several meth- 
ods can be considered. First, utilizing three-dimensional curved channels may lead to elastic instability 
because the pressure hardly suppress spanwise secondary flows for the lack of the degree of freedom. 
Spanwise secondary flows, for example Gortler vortices, are actually observed in three-dimensional 
curved- wall flows |43| . A second normal stress difference in three-dimensional systems also influences 
flow behavior of curvilinear flows of viscoelastic fluids [42] though it can work as both stabilizing and 
destabilizing effects [31| . Second, subcritical instabilities may occur against well-selected finite-amplitude 
perturbations in wavy- walled channels as in the earlier studies of plane- Poiseuille flows [351 US] . The 
well-selected initial conditions are generally produced using the most unstable modes derived in linear 
stability analyses [351 ■ The initial conditions probably have the self-enhancing characteristic that the 
elastic effects store energy in high shear regions and dump it in less sheard regions in order to sustain 
perturbations [45] , 

Finally, we summarize our study. We have numerically investigated laminar steady flows of viscoelas- 
tic fluids in canonical wavy-walled channels at moderate Re and Wi with spectral element method. The 
steady solutions of laminar flows are classified into three groups. The three groups correspond to three 
regions in the Re-Wi parameter space. The steady solutions in the three regions have remarkably dif- 
ferent flow patterns characterized by jet flows near the walls and the vortex detached from a bulk flow 
called separation vortex. The separation vortex and the jet flow near the bottom wall emerge at low Wi. 
The vortex is temporarily vanished or becomes very small at moderate Wi, while the vortex reappears 
at large Wi. The jet flow near the bottom wall disappears and another jet flow appears near the upper 
wall at moderate Wi, while the jet flow near the upper wall gradually reaches to the bottom wall with 
Wi. The separation vortex has been observed in numerical simulations and experiments of both New- 
tonian and non-Newtonian fluid flows in wavy- walled channels [201 SI] ■ The results we have obtained 
are consistent with those of the earlier studies in the point that both the convective and elastic terms 
influence on the flow pattern formation |47| . 

When the flow pattern changes, the viscous wall friction drastically decreases. The drag reduction of 
laminar flows also occurs in experiments and numerical simulations of viscoelastic fluids in curved pipes 
[42l |48] . The flow pattern transition also accompanies the sudden growth of the elastic wall friction. 
The sudden growth of the elastic wall friction is caused by the influence of the elastic stress augmented 
in the bulk flow on the walls. The growth mechanism of the elastic wall friction has similarity with that 



28 



(a) -1.4x10" 



-1.5x10" 



-1.6x10" 



■5 -1.7x10" 



-1.8x10" 



-1.9x10" 



-2.0x10" 



(c) 




10 15 20 

normalized time T 

(c) -1.85x10"' 
-1.86x10"' 
-1.87x10"' 



-1.88x10"^ 



-1.89x10"' 



-1.90x10"' 



-1.91x10"' 



-1.92x10"' 



30 




1.849x10" 



-1.882x10" 



(b) 4.25x10"' 



4.20x10"" 



4.15x10"' 



4.10x10"^ 



4.05x10"^ 



4.00x10"' 




Wi=8.0 




10 15 20 

normalized time T 



12 3 4 5 



10 



15 



20 



25 



30 



normalized time T 



Figure 23: (Color online) Relaxation processes of the velocity component Ux against small random 
perturbations in (a) the convective region Wi = 2.0, (b) the transition region Wi = 5.2, and (c) the 
elastic region Wi = 8.0. The inset of (c) shows the relaxation process magnified in the initial short 
period at Wi = 8.0. Normalized time T is defined as T = t/Wi at each Wi. Horizontal lines denote the 
velocity component Ux of the initial conditions. 
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of the wall friction in Newtonian turbulent flows in which the Reynolds stress t^ = UiUj intensified near 
walls influences on the walls where iti is a fluctuation component of velocity it = u + u. The similarity 
between laminar non-Newtonian fluid flows and Newtonian turbulent flows has been long pointed out 
|49) . The similarity has been also interested in the recent studies in the context of elastic turbulence 

[HIIZI. 

We qualitatively described the mechanism of the flow pattern transition by utilizing the local force 
relation among the centrifugal force, the elastic force, and the pressure gradient on the basis of the 
variation of the terms in the momentum equation. Through the all regions, the streamwise elastic force 
as well as the hoop stress is the origin of the flow pattern formations. When a convective effect and an 
elastic effect balance, the abrupt pattern transition occurs. 

By the scaling analysis, we shown that the stress component is proportional to the Reynolds number 
on the boundary of the first transition in the Re-Wi space. This scaling coincides well with the numerical 
result. The critical Reynolds number Rec ^ 70 obtained in the scaling analysis is reasonable as a 
characteristic value at which an inertial effect becomes considerable |13|. 

We perform the computations of purely elastic flows at Re = 0. Due to the vanishment of the 
centrifugal force, the separation vortex disappears at low Wi. The separation vortex emerges at moderate 
Wi as in the case of the moderate-Re calculations. The vortex vanishment at low Re has been indeed 
observed in the computations of Newtonian fluid flows [5D] . Note that the Weissenberg number at which 
the separation vortex emerges is close to the Weissenberg number at which elastic turbulence starts in the 
experiments of viscoelastic fluid flows for very small Reynolds number Re ^ 1 in curved-wall channels 
[T^ . The modulation of the separation vortex may be related to elastic instability as observed in the 
earlier experiments 1411 . 

We finally perform the computations of Poiseuille-type flows driven by a constant body force. The 
transition of the flow pattern and the growth rate of the elastic wall friction are also observed in the 
Poiseuille-type flows as in the Couette-type flows. 

In the future work, we plan to perform nonlinear stability analyses by utilizing the steady solutions 
which we have obtained in this study in order to investigate a detailed bifurcation diagram in two- 
dimensional wavy-walled channels. In addition, we will employ the three-dimensional direct numerical 
simulations to confirm elastic instabilities in the three-dimensional viscoelastic wavy-walled systems. 
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